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In this paper we present results of two novel experimental methods to investigate the collisional behavior of individual macroscopic 
icy bodies. The experiments reported here were conducted in the microgravity environments of parabolic flights and the Bremen 
drop tower facility. Using a cryogenic parabolic-flight setup, we were able to capture 41 near-central collisions of 1.5-cm-sized 
ice spheres at relative velocities between 6 and 22 cms"'. The analysis of the image sequences provides a uniform distribution of 
coefficients of restitution with a mean value of £ = 0.45 and values ranging from e = 0.06 to 0.84. Additionally, we designed a 
prototype drop tower experiment for collisions within an ensemble of up to one hundred cm-sized projectiles and performed the 
first experiments with solid glass beads. We were able to statistically analyze the development of the kinetic energy of the entire 
system, which can be well explained by assuming a granular 'fluid' following Haff's law with a constant coefficient of restitution of 
£ = 0.64. We could also show that the setup is suitable for studying collisions at velocities of < 5 mms"' appropriate for collisions 
between particles in Saturn's dense main rings. 
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Saturn's rings are one of the most tremendous astronomical 
objects in our solar system, and have been studied for centuries 
by ground-based astronomy (G. Galilei, J. D. Cassini). More 
recently, deep-space missions (e.g. Voyager and Cassini) have 
allowed for in-depth investigations of these structures. Today it 
is clear that the dense main rings (A, B and C ring and Cassini 
Division) consist of myriads of particles orbiting Saturn in Ke- 



and Hameen-Anttilaj ( |1978| l, based on earlier work by |Trulsen| 
( [1971, |1972) . These theories, which were applicable to thin 
granular gases of a single species (i.e. one particle size), were 



extended to dense granular disks ( 


Borderies et al. 1 985 1 Araki 


and Tremaine 1986 Araki 1991 


1, and subsequently explored 


with N-body simulations (e.g. Wisdom and Tremainej |1988[ 



Salo" '1991a"b^ that also included the particle size distribution 



described by,Zebker et al. ( 1985 1. In addition, visco-elastic col 



plerian orbits creating a disk only a few meters thick ( Tiscareno 



2008 1 and radio occultation data (Tyler et al. 1981 Marouf 



etal. 



^FaLl |2007| [Hedman et al] |2007| [Porco et al.| |2008|l. Spec 



lision models, developed by Spahn et al. ( 1995 1 and Brilliantov 



et al. ( 1996 1, have been extended to also treat fragmentation and 
coagulation ("Spahn et al.','2004',iAlbers and Spahn] [2006|. The 



troscopic observations (|Poulet et al.| [2003] |Nicholson et al.| ^^^^^^ ^^^^-^^ (2ooi ) and Salo et ~ 



(2001ishowfliat 



1983 Zebker et al. 1985 1 have confirmed that the rings 



consist predominantly of bodies of almost pure water ice, with 
sizes between 1 cm and 10 m. Many of the observed ring fea- 
tures are caused by interactions with nearby moons and moon- 
lets which gravitationally perturb their trajectories and thus in- 
crease their orbital eccentricity and the ring thickness. Coun- 
teracting this, frequent inelastic collisions at very low rela tive 
velocities (v < 0.5 cm s ', see review by Esposito 2002 1 are 



means to dissipate kinetic energy and thereby re-circularize the 
orbits and confine the height of the ring plane. 

Several theoretical and numerical studies have been carried 
out to investigate how structures inside the dense rings are cre- 
ated, how ring-particle orbits evolve over time and the expected 
outcome of typical collisions between ring particles. Kinetic 
theories were first developed by Goldreich and Tremaine ( 1978[ l 
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the amount of energy dissipated through collisions determines 
the stability of homogeneous dense rings, and thus is a crucial 
parameter in understanding the formation of structures, insta- 
bilities and wakes as observed in Saturn's rings (Colwell et al.[ 
|2006] [20071 [Hedman et allpOOT) [Thomson et al.t|2007j ). 

In most N-body simulations and kinetic theories, the inelas- 
tic collisional behavior is described by a single parameter, the 
coefficient of restitution s, which is the ratio of the relative ve- 
locities v' after the encounter to the relative velocities v before 
the encounter, and is a measure of the dissipation of the ki- 
netic energy. Although ice-paricle collisions play a vital role 
in determining the stability and evolution of Saturn's rings, as 
well as in the formation of icy planetesimals in protoplanetary 
disks, their collision properties have hardly been investigated 
in the low velocity regime. Previous experiments carried out 
by [Bridges et"aL[ ([1984) and [Hatzes et aL[ ( ,1988 ), used a disk 
pendulum to study central impacts of large ice bodies onto a 
flat solid ice target, and showed a decrease in the coefficient 
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of restitution with decreasing impact velocity and increasing 
amount of surface frost. The experiments were extended to in- 



clude grazing collisions of polished ice spheres (Supulver et al 



1995| l. Higa et al. ( 1996 1998, ) conducted impact experiments 
of different-sized ice spheres (r = 0. 14-3.6 cm) into a larger ice 
target. They found that the coefficient of restitution decreases 
with increasing particle size and is independent of temperature. 



A major limitation of all these experimental studies (Bridges 



|etaLl[T984l|Hatzes et al.| [T9881 |Higa et al.|[T996l |1998) is that 
they only treat impacts into a large body of infinite mass and 
infinitesimal surface curvature. In the case of the pendulum ex- 
periments an additional problem is that the mass and moment 
of inertia of the ice sphere are not realistic as they reflect those 
of the sample-pendulum-system instead of that of the individ- 
ual projectile. However, the results obtained by [Bridges et al.| 
( |1984j are commonly used in numerical simulations of Saturn's 
rings. 

In this work we present two novel experimental methods for 
investigating collisions of Saturnian ring particles. In Sect. |2] 
the two experimental setups for collisions of individual ice par- 
ticles and ensembles of particles are described. The results of 
our microgravity experiment campaigns are presented in Sect. 
l3] and their application to Saturn's dense rings is given in Sect. 

al 

2. Experimental setup 

For our microgravity studies we used two different experi- 
mental approaches. The parabolic flight setup was designed to 
study collisions of pairs of individual ice projectiles, whereas 
the Bremen drop tower facility experiment is a prototype, built 
to study a large number of low-velocity collisions within an en- 
semble of particles. 

Microgravity is required to be able to investigate low- 
velocity pair-collisions of free particles with an arbitrary im- 
pact parameter (i.e. from central to grazing collisions), because 
it is the only way to access all the translational and rotational 
degrees of freedom. Parabolic flights offer typically 22 s of 
weightlessness per parabola with a residual acceleration of a 
few times 0.01 go, where is Earth's gravitational accelera- 
tion. Drop-tower experiments achieve residual accelerations of 
better than 10"^ go for up to 9 s. 

2.1. Parabolic flight experiment 

Measurements were conducted on the German Space 
Agency's (DLR) 12* Parabolic Flight Campaign, using a setup 
consisting of a cryogenic cooled experiment body, placed in- 
side a high-vacuum chamber. The experimental apparatus is 
composed of a massive copper block which acts as a thermal 
reservoir, with a rotating sample repository covered by a cop- 
per shield, which plays both a protective and cryo-cooling role. 
The entire system is initially cooled to roughly 77 K with liquid 
nitrogen (N2) spiraling through an attached pipe system. Up 
to 34 pairs of cm-sized ice particles can be stored in the sam- 
ple repository, which allows us to collide approximately 1 pair 
per parabola. To accelerate the ice particles into the collision 



volume, two diametrically opposed pistons are used, driven by 
two synchronized DC motors, which accelerate the particles to 
maximum velocities between 3 - 17 cms"'. These are posi- 
tioned outside the cryo-vacuum chamber, and enter the vacuum 
by differentially pumped feedthroughs. Since the projectiles are 
driven towards each other from exactly opposite directions, the 
resulting binary collisions are predominantly centralized, with 
normalized impact parameters, bjR, close to zero, where b is 
the center-of-mass distance of the two projectiles perpendicular 
to the velocity vector, and R is sum of the particles' radii. The 
collisions are recorded with a high-speed, high-resolution dig- 
ital recording system, operating at 107 frames per second (fps) 
and equipped with beam-splitter optics to provide full three- 
dimensional collision information (i.e. impact parameter and 
relative velocity). For a more detailed technical description 
please refer to |Salter et aL] ( |2009| ). 

41 image sequences of colliding ice spheres (15 mm diam- 
eter) were recorded during DLR's 12* Parabolic Flight Cam- 
paign in April 2008. The samples were produced by freezing 
purified water inside silicone molds in a regular kitchen freezer, 
to produce 15mm-sized, hexagonal-ice spheres. After remov- 
ing the molds from the freezer, the spheres were extracted and 
stored inside a bath of Uquid N2 from where they were individ- 
ually loaded into the sample repository, which had been pre- 
cooled, and cooled further to temperatures of ~ 80 K. This 
procedure is required to minimize the amount of frost on the 
samples' surfaces. Although the maximum exposure time of 
the samples to the humid air was less than one minute, some ex- 
periments show tiny amounts of surface material (which might 
be frost) chipped off during the collisions. 

Since liquid N2 is prohibited on board the parabolic-flight air- 
craft, the cooling procedure had to be stopped before take-off. 
Consequently, the thermal block warmed slightly prior to take- 
off, and subsequently during flight, resulting in sample temper- 
atures ranging from 1 30 - 1 80 K. 

2.2. Drop tower experiment 

A quasi-two-dimensional rectangular glass box of 150x 150x 
15 mm-' volume (see Fig. was constructed to investigate a 
large number of collisions within an ensemble of ~ 100 (cm- 
sized) particles. The test-chamber's walls are made of 2 cm 
thick glass to ensure highly elastic collisions between the sam- 
ple particles and the walls. Two sets of up to 50 projectiles, 
stored opposite each other, were accelerated into the test cham- 
ber with an initial velocity of ~ lOcms ' by two motor-driven 
glass bars that afterwards seal these entry holes (Fig. [TJ. During 
the 9 s of microgravit}{^ the particle positions were recorded by 
a high-speed, high-resolution camera operating at 1 15 fps, and 
two overview cameras (25 fps) positioned at an angle of 16.25° 
relative to the chamber, to enable three-dimensional trajectory 
reconstruction. 

Here we focus on the results of one out of seven micrograv- 
ity experiments conducted at the Bremen drop tower in April 



'The use of the Bremen drop tower's catapuh facility almost doubled the 
duration of microgravity compared to a regular drop without using the catapult. 
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Figure 1: Single image of an ensemble of cm-sized glass spheres colliding in 
microgravity. The black dots on the pailicles are marks for later use in the 
determination of particle rotation (see additional online material). 

2008. This was a prototype experiment for future investiga- 
tions that will be equipped with cryogenic facilities and there- 
fore use ice particles. However in this case, the experiment was 
performed with 92 spherical glass samples of 10 mm diameter, 
with a slightly roughened surface. At the beginning of the mi- 
crogravity period, two sets of 32 particles were injected into the 
test chamber where 28 samples were at rest. The experiment 
chamber was not evacuated, but the influence of the air drag 
acting on a sample particle does not significantly affect its ve- 
locity. A smooth glass sphere (10 mm diameter) has a Reynolds 
number Re ^ 1.6 at a velocity of v = 4mms"', and Re x 62 at 
a velocity of v = 10 cms"'. From this the relative velocity loss 
Av/v due to air drag between two consecutive collisions can be 
estimated to be less than a few times lO"-'. 

Standard methods of image processing and particle tracking 
were used to obtain the particles' positions from the recorded 
image sequence. A statistical analysis of all particles' veloci- 
ties as a function of time was performed, which means that - 
although generally possible - we did not treat individual colli- 
sions but concentrated on the velocity development of the entire 
ensemble. In future experiments, the individual colUsions, their 
impact parameters and the excitation of rotational motion will 
be investigated. The results of our analyses can be found in 
Sect. [32] 

3. Results 

3.1. Parabolic flight experiment 

The binary collision experiments performed during DLR's 
12* Parabolic Flight Campaign occurred in free space and so 



Figure 2: Image sequence of two 15-mm-sized ice spheres colliding at a relative 
velocity of 14 cm s"' . The temporal separation between two consecutive frames 
is -j^ s. The images were captured using a beam-splitter optics. Therefore, 
each image frame contains two views of the set of particles, which are separated 
by an angle of 48.8°. 

were truly three-dimensional events. Thus, data from two pro- 
jections of each event captured using beam-splitter optics (at an 
angular separation of 48.8°) were used to determine the particle 
motion (see Fig. [2|l. The two sets of two-dimensional coordi- 
nates were combined with a transformation algorithm to pro- 
duce the sample coordinates in three dimensions, and the tra- 
jectories were determined by linear fits to the data. From this, 
the relative velocities before and after the collision, v and v' re- 
spectively, and the normalized impact parameter, b/R, were cal- 
culated. Additionally, the velocity components normal and tan- 
gential to the colliding surfaces (vj_ and vy, respectively) could 
be extracted. 

From Fig. [3] it is clear that the impact parameters are re- 
stricted to near-central collisions (due to the limitations of our 
experimental setup), ranging from b/R - to 0.5 with a mean 
value b/R = 0.16. However, simple geometrical considera- 
tions show that in reality glancing collisions are the most fre- 
quent ones. Thus the captured collisions with impact parame- 
ters smaller than b/R = 0.5 cover only a fraction of P(i.p. < 
0.5) = (b/R)^ = 0.25 of the statistically occurring impacts in 
Saturn's rings. From the observed impact parameters and co- 
eflicients of restitution no evidence for a correlation between 
these parameters could be found by linear correlation analysis, 
in which we computed the linear correlation coefficient: 

rb/R = ^ _^ , (1) 
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Figure 3: Obtained coefficients of restitution for ice particle collisions at near- 
central impact parameters in the range b/R = - 0.5. There is no evidence of 
a correlation between the two variables. 
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Figure 4: Obtained coefficients of restitution for ice particle collisions in the ve- 
locity range v = 6-22 cm s"' . A linear correlation analysis shows no evidence 
for a correlation between these variables. 



where r^^^ - 0.10 gives the fraction of the data which can be 
explained by a linear dependence. Based on a similar analysis 
of the coefficient of restitution s and the relative impact velocity 
V (r^ = 0.11; see Fig. ^ a correlation of those can also be ruled 
out. In a detailed analysis of the coefficient of restitution in 
normal {sj_) and tangential direction (ey) no evidence for a cor- 
relation of these values with the impact parameter or the impact 
velocity could be found. Elaborate investigation showed that 
also a dependence of the coefficient of restitution on a combi- 
nation of the collision velocity and the impact parameter could 
be ruled out. 

Hence, we can treat the coefficient of restitution data in a sta- 
tistical way. One way to do this, is to plot the cumulative num- 
ber of collisions with a coefficient of restitution < s (Fig. |5^). 
It turns out that this can be fitted by a uniform distribution (see 
dashed curve in Fig. |5^). The fit provides a mean coefficient of 



Table 1: Summary of the experimental results of the coefficient of restitu- 
tion of low-velocity ice-particle collisions. The velocity range was v± = 
4 ... 22 cm ' and I'y = 0.5 ... 8 cm ' , respectively. 







e± 






=11 




A£|| 


0.41 


0.04 0.24 


0.06 


1.08 


0.19 1.11 


0.38 



£^7 £||: mean normal / tangential coefficient of restitution 

(T^, cr^: error of mean normal / tangential coefficient of restitution 

, fTe,, : standai'd deviation of individual measurements of normal / tangential coefficient 
of restitution from mean value 

A£±, Asf. measurement unceitainty of individual measurement of normal / tangential 
coefficient of restitution 



restitution of e = 0.45 and an overall range of e = 0.06 - 0.84 
(dotted lines in Fig. |5^). 

The obtained normal coefficient of restitution ej_ was treated 
in a similar way. Fig. ^ shows that values of Bj^ fall in the inter- 
val . . . 0.82 (dashed lines). The analysis of our dataset yields a 
mean value = 0.41 with an error of the mean cr^ - 0.04 and 
individual measurement errors of Aej_ - 0.06. The standard 
deviation of the individual measurements from the mean value 
was calculated to cr^^ - 0.24 indicating a significant scatter of 
the values for around eT which is not caused by the mea- 
surement uncertainties. We therefore recommend to include a 
flat statistical distribution of coefficients of restitution (see in- 
sert in Fig. [Sj)) rather than a single value to the models investi- 
gating Saturn's rings. However, it is likely that the wide range 
of coefficients of restitution observed in our parabolic flight ex- 
periments reflects the anisotropic nature of the ice samples' sur- 
face properties, especially the surface roughness or previously 
slightly melted areas altering the collisional properties. 

For the analysis of the tangential component ey only exper- 
iments with tangential velocities vy > 5mms"' (36 out of 41 
experiments) were taken into account to avoid ey - v^/vy ap- 
proaching unrealistically large values when the denominator vy 
approaches zero. The mean value of the tangential coefficient 
of restitution was measured to be ey" = 1.08 with an error of 



0.19 and individual standard deviation of cr. 



fill 



1.11, 

which is in fair agreement with the value found by |Supul-| 
ver et al. (1995i. The large individual measurement errors of 
Asy = 0.38 (due to the small absolute values of v^) do not al- 
low any secure interpretation, however. Table[T]summarizes our 
experimental results for the normal and tangential coefficient of 
restitution of binary collisions between cm-sized ice spheres. 

3.2. Drop tower experiment 

To yield the positions of all particles in all image frames, 
the images were convolved with an image of a single sphere as 
kernel. The probabiUty maxima, which determine the particle 
positions, could then be detected automatically. The tracking 
of all particles over time allowed the calculation of the individ- 
ual particle's mean velocitjj^between consecutive images. The 
resulting particle velocities as a function of time show that the 



^The reader should note that although the term 'velocity' generally de- 
scribes a vectoiial quantity we use it for the absolute value of the velocity. In 
all cases treating the velocity as a vector this is clearly stated. 
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Figure 5: (a) Normalized cumulative number of ice particle collisions with co- 
efficients of restitution < £ for all captured collisions in the parabolic-flight 
experiments (crosses). The distribution of £ can be fitted by a uniform distribu- 
tion which is given by the dashed curve. The mean coefficient of restitution is 
at £ = 0.45 (denoted by the dash-dotted line) with an overall range of e = 0.06 
to 0.84 (dotted lines). The diamonds represent the data obtained from binary 
collisions of glass beads in the laboratory mini drop tower The coefficients 
of restitution span a range of £ a; 0.35 to 0.95 (solid curve), (b) Normalized 
cumulative number of collisions with normal coefficients of restitution < Ex 
for all 41 captured experiments. The normal component has a mean value of 
fiX = 0.4 1 ± 0.04 with a standard deviation of individual measurements of 0.24. 
The measurement eiTor of an individual measurement is 0.06 (indicated by the 
error bars). The insert shows a schematic distribution of coefficients of restitu- 
tion that could be included into numerical simulations. 



system equilibrates very quickly within the first 1.5 s. After 
this injection phase, we observed dissipation of kinetic energy 
through inelastic collisions and transformation to rotational mo- 
tion which was visualized by black marks on the surface of the 
glass beads (see Fig. [T]l. 

From the analysis of the particle velocities within the last 
1 .5 s of experiment time, we obtained a median particle veloc- 
ity of 3.5 mm s"' and determined that 50 % of the particles are 
moving with velocities between 2 and 5 mm s"' (see dashed and 
dotted lines in Fig. |6]), which corresponds fairly well to the de- 
sired velocities for studying collisions in Saturn's dense rings 
(Esposito|[2002l l. 

The decay of the mean particle velocity of the entire ensem- 
ble during the experiment duration is displayed in Fig. |7^. The 
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Figure 6: Normalized cumulative number of glass particles with velocity < v 
in the time interval of 7 to 8.5 s of microgravity duration for the drop-tower ex- 
periments. The dashed line denotes the median particle velocity and the dotted 
lines indicate the particle velocities covering 50% of all particles. 



individual measurements were logarithmically binned and the 
mean velocity was computed. The error bars denote the stan- 
dard deviation of the velocity and experiment time measure- 
ments, respectively. Assuming a constant coefficient of restitu- 
tion e(v) = const., the simple kinetic theory for granular flu- 
ids presented by |Haff| p983) indicates that the temporal rms- 
velocity evolution can be described by the following function: 



v(f) = 



1 



-l-(l - e) ■ n ■ cr ■ t 



(2) 



where vq is the initial injection velocity of the ensemble, n is the 
number density and cr - 4nr^ is the collisional cross-section. 
Eq. |2] can be fitted to the binned velocity data, resulting in a 
coefficient of restitution of e = 0.64 (see solid line in Fig. [7^). 
Thereby, the data represent an average over all impact param- 
eters (i.e. central to grazing impacts) and, thus, collisions are 
only analyzed statistically, but not as individual binary encoun- 
ters. The residual of the fit to the data (Fig. |7j)) shows a good 
agreement of both in the interval from 2 to 7 seconds exper- 
iment duration, meaning that the collisional behavior follows 
Haff's law with e = 0.64. The strong deviation within the first 
2 seconds is due to the equilibration phase of the ensemble af- 
ter the sample injection, whereas the clear change in slope after 
~ 7 seconds is a real effect that can either be explained by an 
increased coefficient of restitution towards low relative veloc- 
ities or the onset of clustering (Goldhirsch and Zanettij |1993[ 
Brito and Ernst 1998 , Miller and Luding, 2004j ) . However, the 
investigation of the latter exceeds the scope of this work. 

We performed a detailed analysis of individual binary colli- 
sions among the constituents of the ensemble in the drop tower 
experiment as well as of binary collisions of identical glass 
spheres conducted in a recently built laboratory mini drop tower 
of 1.5 m height. In the latter experiments we observed 12 cen- 
tral collisions in the velocity range 1 - 4cm s '. The mea- 
sured coefficients of restitution span a range from s ^ 0.35 to 
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Figure 7: (a) Velocity decrease of tlie glass particles in the drop-tower exper- 
iment due to the dissipation of kinetic energy as a function of experiment du- 
ration. The data was binned in logarithmic intervals. It is clearly visible that 
after a short period of equilibration the dissipation is in agreement with Haff 's 
cooling law for granular gases jHaff| |1983 i between 2 and ~ 7 s experiment 
duration. The solid line denotes the relation given in Eq. |2]with a coefficient of 
restitution of £ = 0.64. (b) The residual between Haff 's law and the data shows 
a good agreement from 2 to 7 s experiment time. The strong deviation within 
the first 2 s is due to the equilibration phase after the particle injection, whereas 
the divergence after 7 s can only be explained by a change of the coefficient of 
restitution towards low relative velocities, (c) The scatter of the individual data 
points around the fitted curve shows a broad velocity distribution. 



£ ^ 0.95, are randomly distributed (see solid curve in Fig. |5]l, 
and show no correlation with impact velocity. 

In the drop tower experiments with the particle ensemble we 
have no information about the three-dimensional motion of the 
particles, so that we are forced to restrict ourselves on the mea- 
sured two-dimensional projection of the velocity vector before 
and after the collisions. Splitting the velocity vectors into two 
one-dimensional components, we can describe the distribution 
of the absolute velocity values by a one-dimensional Maxwell- 
Boltzmann distribution function (see right curve in Fig. |8] Fig. 
|9]l with a fixed coefficient of restitution of e = 0.77. Using a 
range of coefficients of restitution, i.e. s - 0.2 . . . 0.9, results in 
a very similar velocity distribution function (left curve in Fig. 
|8]l. The shaded areas in Fig. |8]are the result of a Monte Carlo 
simulation of a one-dimensional Maxwell-Boltzmann velocity 
distribution and denote the 2cr-range. 
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Figure 8: The coefficients of restitution obtained by analysis of the binary colli- 
sions of cm-sized glass spheres (crosses) can be fitted by a Maxwell-Boltzmann 
distribution with a constant coefficient of restitution of e = 0.77 (right shaded 
area). Using a range of coefficients of restitution (e = 0.2 ... 1.0) results in a 
very similar velocity distribution function (left shaded area). The shaded areas 
denote the 2iT-standard-deviation derived by a Monte Carlo simulation. 
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Figure 9: The obtained absolute values of the one-dimensional velocities of the 
binary collisions (histogram with lo" eiTor bars) can be fitted by the positive 
branch of a Maxwell-Boltzmann distribution (solid curve). 



4. Conclusion 

In this work, we present first results of two novel experi- 
mental methods to study the collision properties of individual 
macroscopic ice bodies, allowing free encounters without limi- 
tations present in previous experiments. 

The parabolic-flight setup is suitable for pair-collisions at rel- 
ative velocities > 6cms ' and near-central impact angles. For 
normalized impact parameters of b/R - - 0.5 we could not 
observe a correlation between the coefficient of restitution and 
the impact parameter nor with the relative impact velocity of 
ice particles. The fit of a uniform distribution of the coefficient 
of restitution resulted in a mean value of e = 0.45, and a range 
from s - 0.06 to 0.84. The experiments with equal-sized and 
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equal-shaped glass particles with collision velocities of a few 
cms"' and normal impacts also showed a wide distribution of 
coefficients of restitution between e = 0.35 ande = 0.95. While 
for the ice particles one could argue that surface frost might be 
responsible for the scatter in the coefficient of restitution this 
clearly plays no role for the glass beads. Thus, we expect not 
a single value for the coefficient of restitution for the collisions 
in Saturn's rings, but a more or less wide distribution around a 
mean value. Thus, also very elastic and very inelastic collisions 
concurrently take place in Saturn's rings. The experiments per- 

|T988l l 



formed by Bridges et al. ( 1984 1 and Hatzes et al. 



al- 
though investigating central collisions at lower velocities than 
we did - showed only a weak velocity dependence of s above 
1 cms ' with values of s ^ 0.2 and s » 0.5, respectively, and 
thereby enclose the results obtained from our parabolic flight 
experiments. The experiments by |Dilley and Crawford] ( [T996| l 
show a decrease of the coefficient of restitution with decreasing 
projectile mass and with increasing impact velocities. These re- 
sults do not match those obtained from our parabolic flight ex- 



periments. Recent works by GrasselU et al. ( 2009 1 and Sorace 



|etal.| ( [2009 ^) presenting the investigation of granular cooling and 
low-velocity binary collisions of mm-sized spherical particles, 
respectively, report a decrease of the coefficients of restitution 
with decreasing impact velocities. This finding is supported by 
the collision model of |Brilliantov et aL|(2007) . This means that 
it is unlikely that the deviation of our data points from Haff''s 
law (see Fig. [7^) for f > 7 s is due to an increase in the coeffi- 
cient of restitution. Clustering seems to be a plausible explana- 
tion, but needs further exploration. 

Recent N-body simulations by iPorco et al. (2008i were car- 
ried out to create a model of Saturn's A and C ring. This model 
was used as input for numerical light-scattering studies trying 
to match the photometric properties observed by Cassini. The 
authors found that to reproduce their observations with their nu- 
merical codes, they required a velocity-dependent coefficient of 
restitution of ~ 3.5 times lower than that found by Bridges et al.' 
( ,1984j and 5 times lower than those obtained by Supulver et al., 
(Q995JI. These results could not be reproduced by our experi- 
ments. However, the presence of surface frost or regolith layers 
will alter the collisional properties towards more inelastic be- 
havior A similar trend is known from impact experiments of 
cm-sized spheres into lunar and martian regolith analog mate- 
rial ( |ColweU and Taylor|[T999HColweU, 2003 ) and low-velocity 
collision experiments of high-porosity dust aggregates (see e.g. 



Blum and Miinch 1993 HeiBelmann et al. 2007 [ Langkowski 
etal.,2008F 



In a future microgravity campaign the investigation of ice 
collisions will be extended to impacts onto a solid ice target and 
at arbitrary impact angles using a slightly modified parabolic 
flight setup. 

Using the setup developed as a feasibility study for drop 
tower experiments, we were able to achieve collisions between 
cm-sized glass beads at velocities well below 1 cms"', which 
could be fitted assuming a mean coefficient of restitution of 
s = 0.64. The obtained relative velocities are realistic param- 
eters for simulating the colUsion processes in Saturn's dense 
main rings (Esposito 2002[ l so that future drop tower experi- 



ments using a cryogenic setup and ensembles of water ice sam- 
ples will provide further insights into the collision dynamics of 
icy bodies. These experiments will also incorporate beam split- 
ter optics to provide three-dimensional collision information, 
which proved to be of utmost importance for the analysis. 

The observed collisional properties, like the coefficient of 
restitution s, are crucial parameters for theoretical and numeri- 
cal studies ( |Salo etaL] [2001 ; Schmidt et al.| 2001 ) aiming to ex- 
plain observed structures, instabilities and overstabilities in Sat- 
urn's dense rings. Additionally, the experiments can be utilized 
to determine the mechanism of energy equipartition within an 
ensemble of interacting macroscopic grains. For the explana- 
tion of the extremely confined ring thickness and the means to 
counteract perturbations and excitation by dissipation of kinetic 
energy, more detailed investigations shall be conducted to study 
the influence of regolith-covered surfaces as well as the accre- 
tion and chip-off" of debris in collisions of icy bodies. If ring 
particle collisions are dominated by the particles' surface prop- 
erties and not by their icy cores, a regolith layer will clearly 
reduce the coefficient of restitution considerably ([Heifi elmannl 
etaL] [20071 l^lum and Munch] [T993| [Langkowski et aL| |2008] k 



5. Future Work 

In future experiments using a setup similar to the one used 
for the described drop tower experiment we will collide ensem- 
bles of cm-sized ice particles at cryogenic temperatures. Ad- 
ditionally, low-velocity binary collisions of icy samples will be 
studied in our laboratory mini drop tower Both experiments 
will include stereo optics to gain three-dimensional collision 
information giving insight into the coefficient of restitution as 
a function of the impact parameter and the vectorial collision 
velocity as well as studying the equipartition of translation and 
rotational energy. The sample properties will also be varied 
to investigate the eff'ects of surface properties (frost or regolith 
layers) and shape. 
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Movie caption: High-speed image sequence (115 fps) of an ensemble 
of cm-sized glass spheres colliding in microgravity. The velocity de- 
crease due to inelastic collisions is in agreement with Half's cooling 
law for granular 'fluids' (Half, 1983; ! with a mean coefficient of resti- 
tution of e = 0.64. The black dots on the particles are marks for later 
use in the determination of particle rotation. 
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